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We study the transport coefficients, damping rates and mean free paths of soft fermion coUective 
excitations in a hot fermion-gauge-scalar plasma with the goal of understanding the main physical 
mechanisms that determine transport of chirality in scenarios of non-local electroweak baryogene- 
sis. The focus is on identifying the different transport coefficients for the different branches of soft 
collective excitations of the fermion spectrum. These branches correspond to collective excitations 
with opposite ratios of chirality to helicity and different dispersion relations. By combining results 
from the hard thermal loop (HTL) resummation program with a novel mechanism of fermion damp- 
i-j^ I ing through heavy scalar decay, we obtain a robust description of the different damping rates and 

mean free paths for the soft collective excitations to leading order in HTL and lowest order in the 
Yukawa coupling. The space-time evolution of wave packets of collective excitations unambiguously 
^^ , reveals the respective mean free paths. We find that whereas both the gauge and scalar contribution 

Cn ■ to the damping rates are different for the different branches, the difference of mean free paths for 

both branches is mainly determined by the decay of the heavy scalar into a hard fermion and a 
soft collective excitation. We argue that these mechanisms are robust and are therefore relevant for 
non-local scenarios of baryogenesis either in the Standard Model or extensions thereof. 
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I. INTRODUCTION 
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Q,' One of the fundamental problems confronting particle astrophysics is that of baryogenesis, the origin of the abun- 

I \ dance of matter over antimatter in the observed Universe. The Standard Model (SM) |]1| and extensions thereof offer 

jr~^' the tantahzing possibihty of providing an explanation for baryogenesis with physics at the electroweak scale (for a 

(— I ■ description of mechanisms of baryogenesis see Ref. |g] ) . Although the Standard Model satisfies the three criteria for 

• • . baryogenesis: (i) violation of baryon number, (ii) violation of C and CP, and (iii) departure from thermal equilib- 

. !^ ' riuni &fll , a wealth of theoretical evidence combined with experimental bounds on the CP violating phases in the 

k^ , Kobayashi-Maskawa (KM) matrix, and the current bounds on the Higgs mass, mn > 88 GeV, from LEP2 ||] and the 

;_! ■ most recent lattice simulations M seem to lead to the conclusion that the minimal Standard Model |q| may not be 

5t ! sh\e to explain the observed baryon asymmetry but perhaps extensions could naturally lead to such an explanation 

(for recent reviews see Rcfs. p 10 1). 



An important non-equilibrium ingredient for a successful scenario of electroweak baryogenesis is a strongly first 
order phase transition that proceeds via the nucleation of the true (broken symmetry) vacuum in the background 
of the false (unbroken symmetry) vacuum plus non-trivial topological field configurations, the sphalerons, that are 
responsible for the baryon violating processes H . 

Amongst the several proposals for explaining the baryon asymmetry, we focus on transport aspects related to 



non-local or charge transfer baryogenesis 
electroweak theory (or extensions thereof) 



This scenario assumes a strongly first order phase transition in the 
, pl| . As the nucleated bubbles of the broken phase in the background of 
the unbroken (false vacuum) phase expand, a flux of fermions from the unbroken phase scatter off the bubble walls, 
CP violating interactions near the wall are converted to an asymmetry in the baryon number via sphaleron processes 
in the unbroken phase [g|,|| . 

An important ingredient in this scenario is the transport properties of chiral fermions in the regime of small spatial 
momentum y,g| but in a plasma at temperature T « 100 GeV. In particular there are two different regimes for 
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non-local baryogenesis depending on whether the mean free path of the quarks is smaller or larger than the width of 
the bubble walls 0. 

The hot plasma modifies dramatically the properties of fermions [[L3[-p5[ , in particular when g T 3> Af with g the 
gauge coupling constant and M the zero temperature fermion mass, the fermion spectrum is characterized by two 
branches of collective excitations, one branch corresponds to the in-medium renormalized fermion, with a positive 
ratio of chirality to helicity, and the other branch is a novel collective excitation in the medium, the plasmino, with a 
negative ratio of chirality to helicity ||l^,0. The two branches have a gap Afg// ~ gT corresponding to an effective 
fermion mass that does not break chiral symmetry [[l3|-|l5|. 

The fermion and plasmino branches describe excitations with different chirality properties and therefore these are 
bound to play an important role in charge transfer (non-local) mechanisms of baryogenesis [^|jg] . For soft momenta 
k < gT these branches are very different from each other, while for large momenta k > gT they approach the usual 
fermion dispersion relation, and the amplitude of plasmino excitations (the wave function renormalization) vanishes 
exponentially fast [O-Q . Whereas dispersion relations of soft collective excitations in the plasma within the context 
of baryogenesis had been studied in Ref. [pi and extended to describe neutralinos in the medium in supersymmetric 
extensions in Ref. [|l9| , the damping or thcrmalization rate has only been studied at zero spatial momentum, i.e., the 
collective excitations are at rest with respect to the plasma pQ] or for hard spatial momentum M] k >T. 

Early attempts to estimate the diffusion coefficient made use of Boltzmann equations and a series of approximations 
including a simplification of the gluon propagator [0,n2| . Some of those results were reproduced in pl| for quarks of 
hard momentum , k>T . 

The two limits that were previously studied, zero spatial momentum and very large spatial momentum {k > T), 
do not reveal important features of the collective excitations which are only present at non-zero but soft momenta 
< k < gT where the dispersion relations are different [p"3|-p^ . At zero spatial momentum, both branches of the 
fermionic dispersion relation, i.e., fermions and plasminos coincide and therefore the damping rates for both collective 
excitations must be the same. In the opposite limit k > T, the fermion branch approaches the dispersion relation of 
an ordinary vacuum fermion losing all features of the collective behavior, and the plasmino branch decouples from 
the spectrum because its wave function renormalization vanishes exponentially. 

In this article we present a detailed study of the relaxation rates (damping or thermalization rates) and mean free 
paths for soft collective excitations with non- vanishing momenta k < g T in a theory that includes both gauge and 
Yukawa interactions. To our knowledge such study had not been pursued before within the context of gauge and 
scalar (Yukawa) interactions. Furthermore, since the different collective excitations have different group velocities, 
it is important to understand in detail their mean free paths. This is very relevant for non-local baryogenesis [0,||, 
since the mean free path determines the attenuation length for the transport of chirality by the different collective 
excitations. As emphasized above, the different branches correspond to collective excitations with opposite ratios of 
chirality to helicity, and different dispersion relations. 

Our main observations in this article are that the transport properties of these two different collective excitations 
for soft momenta will in general be different. Since these excitations carry chirality, different damping or relaxation 
rates will result in different mean free paths for collective excitations carrying opposite chirality. As has been argued 
repeatedly in the literature 0Jq,0| the transport of chirality plays a fundamental role in non-local baryogenesis, thus 
a consistent microscopic study of baryogenesis in these scenarios must necessarily address the transport of chirality 
by the soft collective excitations and is the motivation for our study. 

An important aspect of soft excitations in a plasma is that their consistent description requires a non-perturbative 
resummation scheme, the hard thermal loop (HTL) program 123] . For soft excitations with k < gT and vacuum masses 
<^ gT (with g the generic gauge coupling constant) the leading corrections from the thermal self-energies in the high 
temperature limit are of 0{gT), hence of the same order as the tree level terms in the propagators. The leading HTL 
correction is determined by the hard momentum region in the self-energy loop |l4,nq,E3|. However, the damping rate 
of soft collective excitations require that the fermion and gauge bosons propagators and vertices in the self-energy be 
HTL resummed p^ pg] . The HTL resummation rendered infrared finite the damping rate of collective excitations at 
rest in the plasma |23|] , however the damping rate of moving collective excitations is infrared divergent arising from 
the exchange of transverse gauge bosons p4]. Whereas in QCD (and in the SU(2) sector of the electroweak theory) 
the putative non-perturbative magnetic mass provides an infrared cutoff |24| , |25| , in QED (and certainly the U(l) 
sector of the electroweak theory) such a possibility is unavailable and the infrared divergence remains p4| , p5[ . The 
fermion propagator in QED has been studied within the eikonal (Bloch-Nordsieck) approximation in real time P6| . 
This analysis revealed an anomalous relaxation that explained the nature of the infrared divergences in the abelian 
theory and determined that the relaxation time scale of fermionic collective excitations oc [g^T In l/g]~^ ||2^. This 
anomalous relaxation was recently reproduced via different methods in scalar QED |2^ that has the same leading 
HTL structure as QCD and QED for the fermion and gauge boson propagators p8[ . Although there is no magnetic 
screening mass in the abelian theory, the relaxation scales obtained by Refs. [E6|J2W are similar to those obtained in 
QCD Mpl. 



More recently we have obtained the damping rate for fermions in a fermion-scalar plasma and pointed out that 
the decay of a heavy scalar (Higgs) into fermion pair results in a quasiparticle width for the fermion, i.e., the scalar 
decay leads to a fermion damping rate [g^ . In this article we combine the results obtained in QCD and QED for the 
damping rate of the collective excitations with the contribution from the scalars to provide a detailed understanding 
of the relaxation scales and mean free paths for the soft chirality carrying collective excitations in the plasma. 

The goals: In this article we study in detail the damping rates and mean free paths of soft collective excitations 
in a plasma of gauge, fermion and scalar particles as befits the scenarios for electroweak baryogenesis. The focus is to 
provide a quantitative study of the different relaxation scales for the soft collective excitations with opposite ratios of 
chirality to helicity. 

However, rather than studying these aspects in the full Standard Electroweak Model or any particular extension of 
it, we focus on the understanding of the robust and generic features of the damping rates and mean free paths of soft 
collective excitations in an abelian plasma of fermions, gauge bosons and scalars with simple Yukawa couplings. The 
main reasons why the results of this simplification are relevant to the physics of the electroweak theory, are that: 



(i) The structure of hard thermal loops is similar in the abelian and non-abelian theories [[13 16 lq , |23| , |24| with the 
differences only in the gauge group factors and strength of couplings. This is manifestly explicit in similar damping 
rates obtained for QCD and QED pjjM]. A very complete study of the dispersion relations for collective excitations in 
the electroweak plasma has been provided in Ref. jg] . The results of that work clearly show that the HTL contributions 
to the self-energies are universal in the sense that the strong, weak and Higgs sectors yield contributions that are 
similar in form and only differ in the appropriate gauge group structure factors. 

(ii) The simple physics that leads to the damping of fermionic excitations via the decay of a heavy scalar [M is 
fairly robust and transcends a particular model, it is mainly described by the kinematics of the decay of a heavy scalar 
into fermion pair in the medium. 

To model the set of parameters of the electroweak theory we assume a scalar of mass m « 100 GeV according with 
the current LEP2 bounds on the Higgs mass 0, and a first order phase transition temperature T « 100 GeV which 
results in ni w T. This will become important in the analysis of the scalar contribution to the fermion self-energy. 
Furthermore we will focus on the lightest quarks and leptons for which gT ^ M with M being the fermion vacuum 
mass. Only for the lightest fermions we expect soft collective excitations in the medium, since for fermions with 
M ^ gT the in-medium corrections to the self-energy are perturbatively small for weak couplings. For these lightest 
generations, the Yukawa couplings y ^ g and the scalar contributions to the fermionic self-energies are perturbative. 
However, as we will see in detail below, the scalar contribution to the damping rates and mean free paths is important. 

Main Results: The main results can be summarized with the statement that for soft collective excitations that 
carry opposite ratios of chirality to helicity, i.e., fermions and plasminos, the damping rates r-i-(fc) and mean free paths 
X±{k) are different. The difference in these relaxation scales is maximal in the range of momenta < k < gT. The 
difference in the mean free paths of the soft collective modes is dominated by the scalar contribution that we explicitly 



compute in Eq. (3.36). The gauge boson contribution to the damping rate of moving collective excitations arises mainly 
from the absorption and emission of soft transverse gauge bosons for which a hard thermal loop resummation of the 
propagators and vertices is required. On the other hand the scalar contribution to the damping rates arises from 
the decay of the scalar into a hard single particle excitation and a soft collective excitation. Both contributions are 
different for the different branches, in the gauge case the difference is determined by the different group velocities. 

We also provide a wave packet analysis to extract the mean free paths of the moving collective excitations. We 
find that the difference of mean free paths for the collective excitations with different ratios of chirality to helicity 
is mainly determined by the scal ar contribution to the damping rates. The mean free paths and their difference are 



explicitly given by Eqs. (4.9-4.10) 



The article is organized as follows. In Sec. II we introduce the renormalized, real-time Dirac equation for fermions 
in the medium, the renormalization aspects are important for consistency because we need to isolate the wave function 
renormalization factors for the collective excitations. The main reason for studying the effective in-medium Dirac 
equation in real time is that this formulation allows us to extract unambiguously the mean free paths or attenuation 
lengths by focusing on the time evolution of wave packets of collective excitations. This treatment bypasses the 
assumptions of transport by diffusion. Furthermore, in this manner we can make direct contact with the results of 
Refs. [ 2q , p7| , p9| and adapt them to the situation under study. In this section we obtain the different contributions 
to the self-energy and renormalize the ultraviolet divergences. Sec. HI is devoted to a detailed study of the damping 
rates and we establish some of our main results. In Sec. IV we study the time evolution of wave packets of collective 
excitations and extract the mean free paths or attenuation lengths. We then show that the main contribution to 
the difference between the mean free paths for the two soft collective excitations arises from the scalar contribution 
to the damping rate. Sec. V summarizes our results and describes their implications for the standard model (or its 
extensions) and points out the possible modifications to our results in extensions of the Standard Model that include 
light scalars. 



II. EFFECTIVE DIRAC EQUATION IN THE MEDIUM 



An unambiguous identification of the damping rates requires an analysis of the retarded propagators in the complex 
frequency plane. The lifetimes (damping rates) are identified as complex poles in the unphysical Riemann sheet. 

Rather than computing the self-energies in imaginary time and perform the analytic continuations to obtain the 
retarded self-energies in the complex frequency plane, we obtain the effective Dirac equation directly in real time. 
This alternative approach has two main advantages: (i) it allows us to make contact with the real-time results for the 
relaxation of the single particle Green's function presented in Refs. [Eq,p^, (ii) by following the space-time evolution 
of an initially prepared wave packet of collective excitations we can extract directly the mean free paths. 

As mentioned in the introduction we consider an abelian theory with one fermion species and one real scalar coupled 
to the fermion via a Yukawa interaction. We will work in the Coulomb gauge in which W ■ A = because the HTL 
self-energies have a more clear interpretation in this gauge. 

The self-energies in leading order in HTL have been proven to be the same in Coulomb gauge and in covariant and 
axial gauges [ p^p3|j2^ . Ref. |18| provides a discussion of the gauge parameter independence of the HTL contribution 
to the self-energies. 

The Lagrangian density is given by 
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where the Grassmann valued source terms were introduced to obtain the effective Dirac equation in the medium by 
analyzing the linear response to these sources [Q . We have not included the possible gauge couplings of the scalar 
since these will not contribute to leading order in the hard thermal loops and lowest order in the Yukawa couplings. 

The self- interaction of the scalar field accounted for by the term £i[(j>] need not be specified to lowest order. 

We now write the bare fields ^, (j), At, A° and rj in terms of the renormalized quantities (hereafter referred to with 
a subscript r) by introducing the renormalization constants: 
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With the above definitions, £ can be expressed as: 

{l = ^r(i^- ,97"^" + 57 • Kt - 2/</'r - m) *, 

where m{T), M, g and y are the renormalized parameters, and the counter-term Lagrangian Cc.t. is given by 
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and 6£ri[4'r] are the counterterms to be determined consistently in the perturbative expansion by choosing a renor- 
malization prescription. The main reason for going through the renormalization is to extract the wave functions 
corresponding to the collective excitations since these will be important for the calculation of the damping rates. By 
choosing the wave function renormalization Z^ to be defined on-shell at zero temperature, we isolate the finite tem- 
perature from the zero temperature renormalizations unambiguously and recognize the wave function renormalization 
for the collective excitations directly. 



We have included the finite temperature corrections to the scalar mass in the renormalized Lagrangian and explicitly 
introduced a finite temperature counterterm. The scalar mass receives finite temperature corrections also from hard 
thermal loops, proportional to T^. As it will become clear below we are interested in the kinematic decay of the Higgs 
into the fermionic collective modes which will contribute to their mean free paths. By including the finite temperature 
corrected Higgs mass in the Lagrangian allows to treat the kinematic decay in terms of the in-medium corrected mass 
directly. The counterterm 5m{T) will be chosen consistently to cancel the finite temperature corrections to the mass. 
There is a small caveat, however, we envisage our analysis to be valid both in the (homogeneous) unbroken and the 
broken symmetry phases in the case of a strongly first order phase transition. The mass term in the Lagrangian should 
then be defined as the position of the pole of the finite temperature propagator for the Higgs. We will implicitly use 
this definition of the Higgs mass in what follows. Furthermore in our analysis we will assume that in either case, 
unbroken or broken symmetry phases the effective Higgs mass is large enough that its kinematic decay into the soft 
fermionic collective excitations is allowed, in particular we assume that m(T) ^ gT. In the broken symmetry phase 
and for a strong first order phase transition 'm?{T) « to^(O) with m{0) th T k, 100 Gev this assumption seems 
reasonable, in the unbroken phase whether this condition is satisfied depends on the details and strength of the first 
order phase transition. 

If the effective mass in either phase does not allow the kinematical decay of the in- medium Higgs into soft fermionic 
collective modes then there is no Higgs contribution to their mean free path. 

Estimates based on the finite temperature effective potential would seem to suggest that the Higgs mass may not 
be much heavier than gT^ [^0 31 1, but conclusions based on effective potentials in gauge theories are fraught with 
ambiguities such as the gauge dependence of the effective potential which translates into gauge dependence of off-shell 
quantities, and the fact that the effective potential becomes complex [p2|-|34| near its minimum. A Higgs mass defined 
as the second derivative of the effective potential at its minimum is not an on-shell quantity p3] (since the effective 
potential is defined at zero four momentum) , only the position of the pole of the Higgs propagator is an on-shell and 
therefore gauge invariant quantity. Given the uncertainties on the scalar sector of the Standard Model, only a sound 
experimental bound can validate or rule out the assumption used in this work on the value of the Higgs mass. 

We must also point out that our analysis is valid in the unbroken or the broken symmetry phases (under the 
proviso of the kinematic constraint for the decay of the Higgs) but certainly not near the wall of a nucleating bubble. 
The nature of the soft fermionic collective excitations in the presence of a wall is not well understood, since the 
hard thermal loop resummation that leads to the collective excitations is only valid in an homogeneous equilibrium 
background. Thus before we can extrapolate our analysis to the case of collective fermionic excitations near a bubble 
wall a consistent study of the effects of the wall on these collective modes must be pursued. As mentioned above our 
goal is to provide a quantitative analysis of the mean free paths of the soft fermionic collective modes at least in the 
homogeneous phases, with particular attention to the difference of mean free paths of the different collective modes. 

The dynamics of expectation values and correlation functions of the quantum field is obtained by implementing 
the Schwinger-Keldysh closed-time-path formulation of non-equilibrium quantum field theory [p5[-p7[. The main 
ingredient in this formulation is the real time evolution of an initially prepared density matrix and its path integral 
representation. It requires a path integral defined along a closed time path contour. This formulation has been 
described elsewhere within many different contexts and we refer the reader to the literature for details [p5|-p7[ . 

The effective Dirac equation is the equation of motion for the expectation value of the fermion field tplx^t) = 
(\I'.r(x, t)).,, in the presence of the sources fj^rj. Since the details and methods to obtain the effective Dirac equation 
in the medium have been described at length elsewhere [2^, we refer the reader to Refs. |2^,^,gj and simply quote 
the final result: 
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and 5]j:(t — t') is the total self-energy, which to lowest order is the sum of the transverse, longitudinal and Yukawa 
contributions respectively, i.e., Sj:(i — t') = ST(i — t') + S^(t — t') + E-(i — t'). The scalar contribution can be found 
in [29], and for the gauge contribution to the fermion self-energy the main ingredients are the real-time propagators 
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The retarded time evolution must be understood as an initial value problem [E9 36 3^ , which is set up consistently 
by taking the source to be switched on adiabatically from t — — cxd and switched off at f = thus providing the 
following initial conditions to the mean field 

i>^it = 0)=^j:iO), ^^{t<0) = 0. (2.8) 

It is convenient to separate the memory integral in (2^3) from t' ~ — oo to t' = and t' = to t' = t, and for this 
purpose we introduce the kernel crtit — t') defined by the relation 
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Since 'ri^j:{t > 0) = 0, we integrate by parts the memory kernel and finally obtain the equation of motion for i > in 
the following form 
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The equation of motion (2.10) can now be solved by Laplace transform as befits an initial value problem. The Laplace 
transformed equation of motion is given by 

hos - J ■ k - M + S^ (i-fos - 7 • fcj - Sai + (Tj:(0) - saj^{s) V'fc(s) = [ijo + iS^la ~ ^ki^)] V'fc(O) , (2-11) 

where ^^is) and o'j:(s) are the Laplace transforms of V'fe(^) ^^^id iyj:{t) respectively: 
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^g(.s) = / dt e-^* V'fe(i) , '^fe(s) = / dt e-^' aj:{t) . 
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Renormalization proceeds by requiring that the countcrtcrms cancel the divergent parts of the self-energy. 
Specifically, the renormalized self-energy contributions are given by 

'^rfc(O) = '^fc(O) - SM -S^^-k, a^g(s) = a^,,(s) - ijoS^ . (2.12) 

The counterterms will be displayed explicitly below. We thus obtain the fully renormalized effective Dirac equation 
in the medium, which is given by 
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The solution to ( 2.13 ) is therefore given by 
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where we have introduced the full renormalized fermion propagator in terms of the Laplace variable s 
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The retarded propagator is obtained via the analytic continuation s = —iuj + 0"*". Stable excitations correspond to 
isolated poles of S{s = —ito + 0"*",^:) in the physical Riemann sheet in the complex uj plane, whereas resonances 
correspond to complex poles in the second (or higher) Riemann sheet. These resonances are quasiparticles with 
lifetimes determined by the imaginary part of the complex pole and the spectral density for the fermion propagator 
features a Breit-Wigner form with a narrow peak in the weak coupling limit. 



A. The self-energy 



Our main goal is to provide a description of damping and transport of soft fcrmionic collective excitations with 
typical momenta k < gT. In this region of soft momenta, the hard thermal loop contribution from gauge bosons must 
be treated non-perturbatively since it is of order g T. In principle if there are HTL contribution from the scalars, 
it should also be treated non-perturbatively. However we focus our study on the light quarks and leptons for which 
the typical Yukawa couplings y <C g, the HTL contribution from the gauge bosons gives the only non-perturbative 
contribution (at least to lowest order). 

Using the non-equilibrium Green's functions for free fields given in Ref. |g^ and Eqs. (2/7), we find to one loop 
order the following contributions to the fcrmion self-energy: 
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where I]-^*^(i — t') to one loop order 0{y^) are given in Rcf. ||29[. Although the transverse self-energy is available 
in the literature in the imaginary time formulation of finite temperature field theory, their form in real time is not 
readily available. It is given by 
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with bjq = \JM'^ + cp, uJk+q = J {k + q) 
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and distribution functions for the fermion and the gauge boson in the loop. 

These contributions to the self-energy are solely in terms of the bare propagators for the internal lines in the loop. 
For the gauge contribution to the self-energy this is only valid for loop momenta k ^ gT, this is the region of the 
loop integral that gives rise to the HTL non-perturbati ve par t. The region of soft loop momenta k < gT requires the 
HTL dressing of the internal propagators and vertices |^,g^,^ and will be discussed in detail later. 

For the scalar contribution in principle one should consider the different regions of the loop momentum, however, 
since the scalar is taken to represent a Higgs of mass w 90 — 100 GeV ~ T the scalar line does not require HTL 
resummation. However as argued above, the scalar self-energy is perturbatively small as compared to the HTL part 
of the gauge boson self-energy. Its real part will yield a perturbatively small correction to the dispersion relation of 
the collective excitations and will be neglected, however the imaginary part will lead to a contribution to the damping 
rate, which is important since in lowest order HTL the collective excitations are stable. We anticipate now and it will 
be argued more forcefully below, that for a Higgs mass of 0{T) the energy conservation constraint for the imaginary 
part requires that for soft external momentum of the fermion, the internal fermion line must carry a momentum at 
least of 0{T) and therefore in order to obtain the scalar contribution to the damping rate the internal fermion line 
does not need HTL resummation. This is a very important point to which we will come back below and it results in 
that the scalar contribution does not suffer from the infrared sensitivity associated with the gauge boson contribution 
to the self-energy. 

Before we discuss each contribution to the self-energy in detail, we proceed to fix the counterterms to render the 
effective Dirac equation finite to one loop order. 



The ultraviolet-divergent parts of <yj:{0) and d-j:{s) entering in the Dirac equation (2.11) are contained in the zero 
temperature contributions. We use dimensional regularization in three spatial dimensions, introduce a renormalization 
scale /C and find 



^^£(0) = - 






h nmte 



^M = -no (^yI^ + ^ j ^ + finite ' (2.22) 

where the finite parts contain the finite temperature contributions plus zero temperature finite terms. The countert- 
erms can then be chosen in minimal subtraction to be given by 

A convenient choice corresponds to on-shell renormalization at zero temperature, since in this case the residue at the 
poles of the propagators yield directly the wave function renormalizations of the collective excitations. 

III. THE IN-MEDIUM FERMION PROPAGATOR 

After renormalization the zero temperature part of the self-energy is finite and perturbatively small, therefore we 
can neglect it and focus on the finite temperature part. 

As mentioned in the introduction, our goal is to describe the transport of chirality by the soft collective excitations. 
Since the typical size of a non-perturbative gauge field configuration is 0{l/g'^T) [||, the contributions to baryon 
violating processes from quarks and leptons with soft momenta < gT is important S . For these soft fermionic 
excitations, the HTL contribution to the self-energy given by the hard momentum region of the loop integral (i.e., 
k > T) is 0{gT) and must be treated non-perturbatively. 

Furthermore in the region near the electroweak phase transition with T w 100 GeV and g « 0.3 — 0.6 one finds that 
gT ^ M and the vacuum mass term can be neglected for the lightest leptons and quarks. 

Therefore: (i) we neglect the zero temperature contribution, (ii) neglect the (renormalized) fermion mass, (iii) sepa- 
rate the HTL non-perturbative contribution from the hard gauge boson exchange from the perturbative contribution 
of soft gauge boson exchange and scalar exchange. Thus the gauge contribution to the fermion self-energy is separated 
into the hard thermal loop part (5S^jj(s) which is given by the hard loop momentum region fc « T, and a correction 
*S^j:(s) pjl which arises from the soft region of the loop integration k < gT. In this region the internal fermions and 
gauge bosons as well as the vertices must be HTL resummed and a detailed analysis finds that this contribution is of 
order g^T p4| (see below). Whereas the HTL contribution only gives an imaginary part below the light cone (Landau 
damping), the g^T correction gives a contribution to the damping rate [E4| both at rest and for a moving fermion. 
Since the Yukawa coupling for the lightest quarks and leptons is y ~ 10^^^ g and the ratio of the scales that enter 
in the scalar loop is m/T « 1, the scalar contribution will be treated perturbatively, along with that of *I]^j:(s). 

Thus to this order, the renormalized full inverse fermion propagator is 

S-\s,k)^ijos-l-k + 6f:^i:{s) + %^{s) + f:K{s) . (3.1) 

It is convenient to write the full fermion propagator in terms of a non-perturbative and a perturbative part 



S-\s,k)^ 


.5jvks,fc) + 5pi(s,fc) 


'-'atpW' ^) ~ 


= ilos~i ■ k + St.^g{s) 


Sp\s,k)^ 


=%.-(-)+ ^:,(^) • 



(3.2) 

The non-perturbative term will determine the position of the poles, i.e., the dispersion relation of the quasiparticles. 
In the HTL approximation it is known that the poles correspond to stable collective excitations, because the imaginary 
part of the self-energy is non- vanishing only below the light cone (Landau damping). Therefore the damping rate 
for on-shell excitations vanishes to leading order in HTL. The real part of the perturbative contribution will provide 
a small shift to the dispersion relation, whereas the imaginary part will provide the damping rate of the fermionic 
collective excitations on-shell. Because the correction to the dispersion relations from the perturbative part is small, 
we will neglect it but instead focus on the imaginary part since it will determine the damping rate which vanishes to 
leading order in the HTL. 

Since the fermion self-energy is a sum of the contribution of the gauge bosons and that of the scalars, we study 
each contribution separately. 



A. Non-perturbative contribution to the fermion self-energy: hard thermal loops and collective excitations 



We now focus on S^p{s, k), which is given by the bare propagator plus the HTL contribution from the gauge boson 
exchange. The HTL approximation is obtained by considering the region g « T in the loop integrals, i.e., the loop 
momentum is hard and the external momentum is soft k < gT. 

The HTL contribution to the fermion self-energy is obtained directly from the Laplace transform of Eqs. ( 2.19| - 2.21 ) 
in the HTL limit. In terms of the Laplace variable s, it is found to be given by 



6Kj:{s) 



9^ 
16k 



70 In 



k 



+ 7 • fc 



*-5 , 

2 In 

k 



The analytic continuation of 5Ti^^{s) in the complex s-plane reads 

JE^j(s = -iuj ± 0+) = (5S^(w) = 5T.p,j:{uj) ± i(5E^j(w) , 



where 



<5S„r(tj) 



<5S.c(w) 



16/c 



2t^2 



70 In 



uj + k 



uj ~ k 



+ j-k 



UJ 

2- -In 
k 



u) + k 



U! — k 






jo-i-kj e{e~Lo^) , 



(3.3) 

(3.4) 

(3.5) 
(3.6) 



and k = k/k. 

To leading order in ff (HTL), the full fermion propagator is given by *S~^{ijJ, k) — W70 — 7 • fc + (5E^r(w) which can 
be written in the form p4|jlq ] 



'S{io,k)^- *A+(c^,fc)(7o-7-fc) + *A_(c^,fc)(7o + 7-fc) 



where 



'■A±^(w, fc) =uj^k 



»'n 



i^-d'i'Q 



±1 



«»(l)^5 



In 



UJ + k 



UJ — k 



iTTQ{k' ~ uj') 



(3.7) 

(3.8) 
(3.9) 



Here, Me// — gT/\/8 is the effective fermion thermal mass induced by the gauge coupling. 

The collective excitations correspond to the poles of the fermion propagator: there are two branches for positive 
energy (and two for negative energy) which are the solutions of *A^^(a;,fc) — 0, the Dirac spinors associated with 
these solutions have opposite ratios of chirality to helicity x ^^'^ the residues at the poles are Z±{k), respectively, 
with |l3|-|l8| 



*Ag^{uj,k) = Q=^uj = uj±{k) , Z±{k) = 



^Kf 



[ujl{k)^k^ 



(3.10) 



where the dispersion relations uj±{k) are shown in Fig. 1. The upper branch corresponds to w+(fc) and describes 
collective excitations with ratio of chirality to helicity x = +1; the lower branch corresponds to uj^{k) describing 
collective excitations with ratio of chirality to helicity x — ^^- The upper branch corresponds to the usual fermionic 
excitation, whereas the lower branch has vanishing group velocity (the minimum in the dispersion relation) at fcmin ~ 
0.4 Mgff and describes a new collective excitation in the medium, it has been named the plasmino to emphasize that 
it is a fermionic excitation that only exists as a collective excitation in the plasma |l3-^,Qq|. The vanishing group 
velocity at fcmm is a novel feature of the plasmino collective excitation and it will be important for the damping rate 
and mean free path. 

We now collect some of the relevant properties of these solutions which will be relevant for the interpretation of the 
results (see Ref. Q for details) 






k » Meff ■■ w+(fc) « fc + 



^4^// 



n^ , 2fc 

aj_ ( A; j w K H exp 
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Z+(fc) ', 



2fc2 

Z-{k) 



1-ln 



2fc2 



2fc^ 



gM2 



■ exp 



e// 



2P 



(3.11) 



Two noteworthy aspects of these expressions wiU be important for the interpretation of particular features of the 
damping rates: (i) the dispersion relations are always above the light cone, (ii) the wave function renormalization 
(residue at the pole) Z-{k) vanishes very fast for k > kmin- Thus the collective excitation with x = ^1 only 
contributes in the region of very soft external momentum k < kmin ■ 

B. Perturbative corrections: the damping rates 

Neglecting the fermion vacuum mass term, the general form of the self-energies is dictated by rotational symmetry, 
hence *S^g(s) and ^^As) can be written as 



%i:{s)^~^SJo*f^f\s)+J■k*f:^^\s) 



;=(°h 



t^W, 



J:K{s)^^Jose'y'is) + J■ke)^>is), 



(3.12) 
(3.13) 



where we have used the notation of Ref. [E9[ for E^-(s). The analytic continuation of these contributions in the 
complex s-plane are defined by 



E^(w) = E^(s = -tu ± 0+) = E^ Juj) ± i E^(w) . 

rk^ ' rk^ ' Rk^ ' Ik^ ' 



To order g^ and y^, the full fermion propagator now reads 

1 



S{uj, fc) = - A+(cc;, fc)(7o - 7 • fc) + A_(cc;, /c)(7o + 7 • /c) 



where 



-^ — 1 



A± (tj, k) ^ *A±'(tJ, k) + n±(a>, k) , 

n±(.^, fc) = iTr [(70 T 7 • fc) i^^rkk^) + ^.jk^) 



(3.14) 
(3.15) 



(3.16) 

(3.17) 
(3.18) 



with * A_|_ ( c^,fc) given by Eq. (3.8) and the real and imaginary parts of n±(a;) can be read off from equations 
(p4|)-(p5|). 

The poles of the full fermion propagator determine the excitations in the medium. The position of the poles are 
obtained from the zeros of A_(_ (w, fc) for lo = ujp{k) — iT{k). In the narrow width approximation r(fc) <^ ujp(k), the 
position of the complex poles are determined by the following equation 



*A±i(a;p(fc),fc) +n±,fl(a;p(fc),fc) -^ [Z±i(c^p(fc)) r(fc) + sgn(r(fc)) Il±j{u;pik),k) 
where we have made use of Eqs. ( 3.14| ) and ( 3.1E| ), and 

Z±{iUp{k)) 



= 



(3.19) 



d*Ag\uj,k) 



diO 



.±(fe) 



are the residues at the poles for the collective excitations. 
The real and imaginary parts of the above equation read 
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*Al\ijp{k), k) + n±^R{ujp{k), fc) = , (3.20) 

r(fc) + sgn (r(fc)) z±{ujp{k)) n±j{LUp{k),k) = o . (3.21) 



To leading order in HTL [i.e., neglecting TI±^r in Eq. ( 3.20 )] the real part of the poles are given by the dispersion 
relations of collective excitations 

^p(fc)-|j"+!S f^""' (3.22) 

^^ ' I ±ijJ_(k) tor plasmmo . ^ ' 



The solution for the imaginary part is obtained by replacing u!p{k) with (p. 22) in Il±j{ujp{k), k) to this order 



However, the equation for the imaginary part (3.21) does not have a solution because Il±j{ujp{k),k) > (see 
Eqs. ( 3.261 ), (3.34) and the discussion below it) and Z±{Ljp{k)) > 0. Therefore, there is no complex pole in the 
physical sheet. This is a well-known result: if the imaginary part of the self-energy on shell is positive there is no 
complex pole solution in the physical sheet, the pole has moved off into the unphysical (second or higher) sheet. If 
the imaginary part on shell is negative, there are two complex poles in the physical sheet corresponding to a growing 
and a decaying solution, i.e., an instability. It will become clear below that Il±j{LLip{k),k) > corresponding to a 
complex pole in an unphysical sheet. In this case the spectral density features a Breit-Wigner resonance shape with 
a width given by the damping rates 

r±(fc) == Z±{k) U±j{Lu±{k), k) , (3.23) 

which in general determines an exponential fall-off of the fermion propagator in real time: the amplitude for collective 
excitations with x = ±1 fall off as 6^'"='='^'''-'*, an interpretation borne out by the real time evolution of the initial value 
problem [ p9|J3qj37| | . An important exception to this analysis appears in the Abelian gauge theory because of infrared 
divergences in the damping rate [£6|,^ . These will be analyzed in detail below. 

1. Gauge boson contribution 

We now focus on the contribution to r-i-(fc) from the soft loop momentum region of the gauge boson exchange. 

For soft momenta fc « ^T, the non-HTL contribution to the fcrnrion self-energy from the gauge boson, i.e., *I]^j(cli) 
requires HTL resummed internal propagators and vertices and is subleading by one power of (7, a thorough study 
is presented in [p3 24|. The main ingredients for a consistent computation of the non-HTL contribution are: (i) 
the HTL resummed internal propagators, (ii) the HTL resummed vertices [p3|-p6[. The resummed propagators are 
obtained from the spectral representation in terms of the HTL spectral densities both for fermions and gauge bosons 
(longitudinal and transverse). When all of the lines at the vertex are soft the differential form of the Ward identities 
from the HTL resummed fermion self-energy can be used to obtain the resummed vertex [g3|-g6[ . We refer the reader 
to Refs. [p3|-p6[ for the details. 

The analysis leading to the expression of the damping rates Eq. ( p. 23 ) relies on the existence and smallness of the 



imaginary parts of the self-energy on the mass shell of the collective excitations. However in gauge theories potential 
finite temperature infrared divergences could invalidate these conclusions. The infrared divergences associated with 
the exchange of longitudinal gauge bosons are a result of small angle Coulomb scattering and at zero temperature 
are those of Rutherford scattering. At finite temperature the longitudinal gauge boson (instantaneous Coulomb 
interaction) is screened by the Debye screening mass niD o<: g T which cuts off the infrared and leads to a finite 
contribution to the damping rate from longitudinal gauge bosons. For soft external momentum this contribution has 
been found to be given by p^ , p6[ 

T^^s(^k)^aAT , (3.24) 

with A a constant that can be found nmncrically [Q and in QCD is of 0(1) and a = g^/47r. 

For a fermion excitation at rest (fermions and plasminos coincide for fc = 0), the damping rate has been com- 
puted ||23|-|25| and found to be given by 

r3{k^0)=aCT , (3.25) 



with C again being a numerical constant of 0(1) [£3|,^ . The reason that we do not specify the constants quantitatively 
is because we are interested in the damping rate for moving collective excitations. When the collective excitations are 
at rest there is no difference in their dispersion relations and therefore their transport properties are identical. 
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In QCD (and also in the SU(2) sector of the Standard Model) the potential infrared divergence arising from the 
exchanged transverse gluon propagator is conjectured to be screened by the non-perturbative magnetic screening mass 
rrimag ~ Q^T {g is the QCD coupling constant) and the longitudinal gluon is Debye screened with m/j w g T. Thus 
the infrared divergences are cured by electric and magnetic screening lengths. The detailed analysis of Rcfs. [2^ J2^] 
lead to the result that the damping rate of moving quasiparticles with momenta k ^ g^T in a non-abelian plasma 
are given by (up to an overall constant that depends on the gauge group structure) 

Tiik) ^ ^\v^ik)\ \n^ , (3.26) 

where v±{k) are the group velocities for the fermion and plasmino branches. This expression is not valid for fc = 
where the damping rates of fermion and plasmino at rest do not have the logarithmic behavior in terms of the gauge 



coupling p3 24 1, nor near k = kmin where the group velocity of the plasmino branch vanishes 



We point out that our results for the mean free paths of the collective excitations and their difference will be 



restricted to the domain of validity of Eq. (3.26), furthermore we emphasize that this regime of validity has been 
pointed out in Refs. ||24|- |26| , |38| ). In the region where the plasmino group velocity vanishes there are contributions to 
its damping rate that must be understood in detail, such a task is clearly beyond the scope of this article. 

Since the HTL structure in QED is similar (up to gauge group factors) to that of QCD (and scalar QED), Pisarski |2J] 
suggested that a similar form of the damping rate should be valid in a QED plasma, despite the fact that there is no 
magnetic screening in the Abelian theory. 

In QED the transverse photon propagator is only dynamically screened via Landau damping and the infrared 
divergences remain, possibly to all orders in perturbation theory. These divergences led to questioning the quasiparticle 
interpretation of moving charged excitations ||3^,Q . More recently [^ a detailed study of the fermion propagator in 
the Bloch-Nordsieck (eikonal) approximation in real time revealed that for lod t \v±{k)\ ^ 1 

Sk{t) w exp[-aT|i;±(A:)| t\n{wD t \v±{k)\)] , (3.27) 

with loj:) ~ gT the Debye frequency and v±{k) the group velocity of the fermion and plasmino branches. Although 
this is not an exponential relaxation that would emerge from a Breit-Wigner resonance shape of the spectral density 
as argued above, it does reveal a particular time scale from which a damping rate can be extracted and is given by 

Tl(k)KaT\v±{k)\\n- . (3.28) 

9 

This result has also been recently found in scalar QED with an alternative method based on the renormalization 
group and within the context of an initial value problem as followed here E^. Since SQED, QED and QCD share 
the same HTL structure to lowest order pq,n8| , the analysis via the renormalization group furnishes an independent 
confirmation of the eikonal approach. 

At this stage it is important to describe the physics that leads to the damping rate from the gauge boson contribution, 
which as it will be seen below is rather different from that of the scalar. The constraints from energy momentum 
conservation in the imaginary part of the self-energy can only be satisfied below the light cone where the HTL 
resummed gauge boson propagator has support that arises from the Landau damping cut. The infrared process that 
leads to the non-exponential relaxation is the emission and absorption of soft photons or gluons at almost right angles 
with the moving fermion [p4| . 

Thus we summarize the contribution to the damping rate of the collective excitations by the HTL resummed gauge 
boson exchange: 

Tl{k)KaT\v±{k)\\n- for k^g'^T, (3.29) 

9 

T%{k)7iiaT for fc = 0. (3.30) 

The damping rate for moving fermions is not yet available for the whole range of group velocities and clearly a 



better understanding of the infrared region must be pursued. In particular Eq. (3.2£) is not valid near the plasmino 
minimum when the group velocity vanishes, since subleading contributions become important in this region [p8| . 
Fig. 2 displays r^(fc)/Afe// vs k/M^ff for g = 0.3. 
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2. Scalar contribution 



In Ref. [ p9| the damping rate for a fermion in a scalar plasma was derived. There it was found that the decay of 
the scalar into fermion-antifermion pairs results in a quasiparticle width for the fermion. 

Here we study a different aspect, that is the lifetime of the collective excitations both fermions and plasminos due 
to the process of scalar decay into the collective excitations. This is one of the novel co ntrib utions of this article. 

Writing the Laplace transform of the scalar contribution to the self-energy as in Eq. ( 3.13 ), the scalar contribution 

to n±(aj, k) in Eq. ( 3.1^ ) are given by [see Eq. ( |3.13 ) for the definition of the e- (w)] 



-(0), 



-(1)/ 



W^{oj,k)=uje'^>{Lj)±ke^y{uj) . 



(3.31) 



-W/ 



The full expressions for the £- (w) to one loop are given in Ref. 



I, we only quote here the most relevant features 



of their imaginary parts to clarify our arguments. The imaginary parts of these coefficients are given by 






J"h 



^0), 



— sgn(.)^p^''(|.|) + p^^(-|.|) 



-sgn(u;) p^}\\uj\) 



4"< 



i"i) 



in terms of the following one-loop spectral densities 
(0)/ N 2 /■ d^l 



Pfi-) 



X [6{lU - LUk+q - Wg)(l + Hk+q 

d^q 1 k ■ q 



y 



(27r)3 2uJk+qUJq k"^ 
X [5{lO - LUk+q - 'ig)(l + Tik+q - flq) - 5{ljJ 



- Wfc^ 



-UJq){nk^ 



(3.32) 



nq) + S{U0 ~ UJk+q + LLlq){nk+q + Uq)] 



(3.33) 



where LOk+q = y mn? + [k + q)'^, and Uk+q and Uq are the respective distribution functions for the scalar and the 
fermion in the loop. The terms that contain the 5{lj — tOk+q + oJq) do not give a contribution to the damping rate. 
These arise from the processes Tp ^^ (j) + ip (with ip and denoting the fermion and the scalar, respectively) and result 
in the usual two particle cuts with support for |ijj| > ^/W+(jti + M)^ (here, M and m are the respective masses of 
the fermion and the scalar in the loop). Thus we only consider the terms proportional to Uk+q + Uq in 11 j_ j{uj,k) since 
only the delta functions that multiply these terms will have support on the mass shell of the collective excitations as 
explained below. These terms lead to the following imaginary part of the scalar contribution 



ni,(cc.,fc) = V 



d^q Uk+q + Uq 



LOq =F Sgn(w)fc • q) 5{\uj\ - UJk+q + LOq) 



(27r)3 2ujq2uJk+q 

-\- \U)q ± Sgn(u;)fc • q] 5{\U}\ + UJk+q - Ulq) 



(3.34) 



It is straightforward to show that IIj. j{ijj, fc) > from the energy conservation constraints: \lij\ ± (D^ =p uJk+q = 0, and 
therefore this imaginary part is not associated with an instability but with a true resonance, i.e., a complex pole in 
an unphysical sheet and its real time interpretation corresponds to an exponential fall off of the amplitude of the 
collective excitation. 

The first delta function determines a cut in the region < \uj\ < y/k"^ + {ra ~ M)^ and originates in the physical 
process (f> —> ip -\-ip whereas the second delta function determines a cut in the region < |a;| < k and originates in the 
process (p + ^ —> ^. The first cut describes the decay of the scalar into fermion-antifermion pairs, the second cut for 
( uj^ < k^) is associated with Landau damping. Both delta functions restrict the range of the integration variable q 
(see below). However, since the dispersion relations for the collective excitations are always above the light cone, i.e., 
U!±{k) > k, the contribution from the second delta function vanishes on the mass shell of the collective excitations. 
Therefore only the first term in Eq. (3.34) contributes to the lifetime of the collective excitations. 

It is at this stage that we are in position to formalize the arguments presented in the previous section to justify a 
one-loop computation of the scalar contribution to the self-energy of the collective excitations in terms of free field 
propagators (non HTL rcsummed). 
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The argument hinges upon two important features of the spectrum of the fermionic collective excita tions: (i) The 
dispersion relation is always above the light cone, this feature guarantees that the second cut in ( p. 34 ) with support 
below the light cone will not contribute to the imaginary part evaluated on the mass-shell of the collective excitation, 
(ii) Under the assumption that the scalar mass m^ gT and an external momentum k fn gT and therefore \uj\ ~ gT 
on the mass shell of the collective excitations, the first delta function is satisfied only for tUq « ^Jq^ + m? ^ gT, hence 
the fermion in the loop is the ordinary single particle fermionic excitation, not the collective excitations associated 
with soft moment a. T his is clear from the properties of dispersion relations and residues for the collective excitations 
displayed in Eq. (3.11) since for large momenta only the fermion branch survives (the wave function renormalization 
of the plasmino branch vanishes) and approaches the vacuum fermion dispersion relation. Hence in summary: the 
contribution from scalar exchange to damping rates of the soft collective excitations is obtained from the one loop 
self-energy in which the internal scalar and fermion lines are the free particle propagators, but the external line are 
the fermion or plasmino collective excitations. 

This is a remarkable result, the fermions and plasminos acquire a width through the induced decay of the scalar in 
the medium and the damping rates for the soft collective excitations are different. For large scalar mass m^ gT 
the kinematics dictates that the scalar decays into a soft collective excitation, either fermion or plasmino, and a hard 
fermion. 

In the limit of M <C 5 T, we finally find 



n|,,(c.,A:) 



y 



327rfc2 



92 ('^) 



<ll{'-^) 



dq [uk+q + Uq) [2kq =F sgn(w) (2 \lo\ q + uj^ - Ll>1)] 



(3.35) 



for |a;| > fc, where tJk — Vm'^ + fe^ and 



iU^) 



■UJf 
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From the expression for the total damping rates Eq. ( p.23| ), the scalar contribution to the damping rates for the 
fermion and plas m ino c ollective excitations is found by evaluating the scalar contribution of the imaginary part of the 
self-energy (3.34)-(3.35) on the mass shell of the collective modes. We find the damping rates to be given by 
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i) 



.-/3'^^ 



(J2(w±(fc)) 



.±(fc)) 



(3.36) 



is the dilogarithm function. ^_^ 

Although these expressions are somewhat unwieldy, there are important features that are very revealing in Eq. ( 3.36 ) : 
these are the wa ve fu nction renormalization factors Z±{k). As we discussed in the previous section and is clear from 
the expressions ( 3.11 ) the wave function renormalization for the plasmino branch vanishes very fast for k > kmin ~ 
0.4 Meff- The damping rates ^±{k) are displayed in Fig. 3 vs k/Meff for g — 0.3 and m = T. The upper solid 
line is the damping rate for the fermion branch, and the lower dashed line is the damping rate for the plasmino 
branch, the rapid decrease for the plasmino branch beyond k > M^ff is a result of the vanishing of the wave function 
renormalization. 



IV. MEAN FREE PATHS: A WAVE PACKET INTERPRETATION 

The notion of mean free paths in a plasma in which different collective excitations with different group velocities 
are present is rather subtle and requires a careful treatment. It is here where the real time description in terms of 
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the in-medium Dirac equation is fruitful. We obtain the mean free paths by studying the space-time evolution of an 
initially prepared wave packet. This formulation identifies unambiguously the mean free paths without resorting to a 
diffusion interpretation of transport. 

In particular we are interested in transport of chirality by the collective excitations, which evolve independently in 
time. Since the collective excitations are eigenstates of the effective Dirac equation in the leading HTL order, they 
are given by the solutions of 



''Az\co, fc)(7o - 7 • fc) + *A-\uj, fc)(7o +j-k) Ms = -i^) = 



(4.1) 



Let us consider that the initial expectation value V^i:(0) in Eq. (2.14) is of the form 



V'^(o) « u^ 



-^{k^kaf 



(4.2) 



where the spinors Ui^ are annihilated by 70 T 7 • ^: i-e. , have chirality to helicity ratios x = ± 1 respectively. This initial 

state describes a wave packet of collective modes localized in momentum at ko- In space it describes a wave packet 
localized over a spatial extent R around the origin with expectation value of the momentum ko- To describe an almost 
monochromatic beam of collective excitations we will choose R ^ uj~^ since oj~^ is the typical spatial extent of the 
quasiparticle states (screening cloud). The quasiparticle is described by a spectral function of the Breit-Wigner form 
(although this is not evident for the logarithmic relaxation in QED, it has been argued to be a good approximation 
in this case [ B6|). The time evolution can then be obtained by inverse Laplace transform, the details of which can be 
found in Rcf. ||37|]. In the narrow width approximation r^(/c) <^ '^±{k), we find the following time evolution of the 
initial wave packet 



iP^{x,t) 



cPke'''-^ Z±{k)Ul 



-i^(fc_fe„)2 



o±(k)t g-r'±(fc)t 



(4.3) 



where r^(fc) is the total damping rate, i.e., gauge plus scalar contribution and Z±{k) the residue at the quasiparticle 
pole (wave function renormalization). Since the initial wave packet is strongly peaked in momentum, we can perform 
the integral over momentum by expanding LO±{k) and T*^{k) around k — ko and using the narrow width approximation. 
We find 



il)^{x,t) (X Z±{ko)U^ e" «*w e '±^-°'-e' 



1^1777 -r^(fc„)t J^ko-{S-Vp,± t) 



(4.4) 



where 



2m _ D2 , ■ d'^^±{k) 



X±{t)^x-Vg^±t , R\t)^R^ + 



Vg^± = ko 

Vp,± = k. 



duj±(k) 



dk 

U!±{ko) 



koV±{ko) 



dk^ 



(4.5) 
(4.6) 
(4.7) 



Obviously Vg^± and Vp,± are the group and phase velocities, respectively. We then see that the peak amplitude of the 
wave packet attenuates in space on a distance scale given by the mean free paths 



X±{ko) 



\v±iko)\ 

r*±(fco) 



(4.8) 



where v±{k) are the group velocities for the fermion and plasmino branches and we took the absolute value because 
there is a region of negative group velocity for the plasmino branch (see Fig. 1) before the minimum at fc « 0.4 M^fj. 
Eq. (4^) is the main reason for the wave packet analysis: it unambiguously identifies the mean free paths of the 
collective excitations in terms of the group velocities and damping rates for each branch of collective excitation 
independently. 

Since t he ga uge contribution to the damping rates of collective excitations is proportional to their group velocities 
(see Eq. ( 3.29| )) we then highlight the fact that whereas the gauge contribution to the damping rates of the collective 
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modes are different, their mean free paths are the same. Hence the difference of mean free paths which is the goal of 
our study is solely given (to this order) by the scalar contribution. This is the main result of this study. 



Using the expression for the gauge contribution to the damping rate of moving quasiparticles given by Eq. ( p.29[ ) 
we can write the mean free paths in the following form 



A±(fc) = ^g 



1 
1 o 



\v±{k)\ n 



aT \n- . (4.9) 

9 



This expression reveals clearly that the difference of mean free paths for the fermion and plasmino collective 
excitations is mostly given by the contribution to the damping rate from the scalar sector. This is one of the 
important results in this article that we want to emphasize: although the largest contribution to the damping rates 
for soft collective excitations arises from the gauge boson contribution to the self-energy, since a S> y^, the difference 
of mean free paths is determined mainly by the decay of the heavy scalar. Fig. 4 shows [A+(fc) — A_(fc)]Me// vs 
k/Meff for g = 0.3, y = 10^'^ and rn = T. This figure displays one of the important results of this article. The 
strong peak near the plasmino minimum is a result of the vanishing of the plasmino group velocity and as discussed 
previously the gauge contribution to the damping rate of the plasmino branch is not reliable near this region but is 
so beyond and before this point. 

In particular away from the minimum from the plasmino branch and using j/^ <C a, the difference in mean free 
paths is approximately given by 



A+(fc)-A_(fc)«^ 



ri(fc) TUk) 



\v-m v+{k) 



(4.10) 



Thus the fermion and plasmino soft collective excitations have different mean free paths and the difference is 
determined mostly by the heavy scalar decay into a hard and a soft collective excitation. 

V. SUMMARY, IMPLICATIONS FOR THE STANDARD MODEL AND CONCLUSIONS 

We have focused our attention on the transport properties of soft collective excitations in an abelian gauge, fermion 
and scalar plasma with the goal of understanding chirality transport as relevant for non-local baryogenesis. 

Our main observation is that the transport properties of soft collective excitations are different for the different 
branches of the fermionic dispersion relation. Since these collective excitations carry ratios of chirality to helicity 
X = ±1, different transport coefficients, i.e., damping rates and mean free paths for these excitations will result in a 
differential transport of chirality. There are two main physical mechanisms that lead to different transport properties: 
(i) The absorption and emission of soft gauge bosons that are dynamically screened by Landau damping on scales 
g T and by a magnetic mass for the non-abelian sector on scales g^T . (ii) The decay of the heavy scalar (Higgs) 
particle into hard fermion and a soft fermion or plasmino results in a contribution to the damping rate of on-shell 
collective excitations. The physical aspects of these two phenomena are fairly robust, in particular the hard thermal 
loop structure of the fermionic self-energy has the same form in the abelian (QED) and non-abelian (SM) theories, 
with the only differences being in the overall scales (see for example Ref. ^ wherein the contributions of the weak and 
scalar sectors to the leading HTL fermion self-energy had been computed in detail). The decay of a heavy Higgs into 
a hard fermion and a soft collective excitation is a consequence of simple kinematics in the heat bath and is therefore 
fairly independent of the model. Obviously the strength of the Yukawa couplings and the group representation of the 
scalar fields will change the overall normalization of the scalar contribution. 

Therefore we believe that the results obtained in this article are relevant for non-local baryogenesis in the Standard 
Model or extensions thereof and that the main ingredients that determine the different transport coefficients described 
in this study are fairly robust. 

We have combined detailed studies of the damping rates of moving fermionic quasiparticles ||23|-|26| and extended 
the recently studied phenomenon of fermionic damping rates via heavy scalar decay |29| to provide a description of 
the damping rates and mean free paths of collective excitations with momenta k ^g^T to leading order in Hard 
Thermal Loop resummation and lowest order in Yukawa coupling. This analysis is valid for the lightest quarks and 
leptons in the Standard Model (certainly not the top quark) and both in the symmetric and broken symmetry phases, 
since for the lightest fermions gT ^ M with M the vacuum masses. The main uncertainties in our results stem from 
the infrared sensitivity of the gauge contributions to the fermionic damping rate, in particular the magnetic sector 
and its non-perturbative screening scale. This is still an ongoing subject of active study ||24|^6| and certainly beyond 



the scope of this article. However as argued persuasively in P^-M the results given by Eqs. (3.28) are trustworthy 
for k ^ g^T and away from the minimum in the plasmino branch when the plasmino group velocity vanishes. For 
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a first order phase transition the Higgs mass does not vary much near the phase transition from its vacuum value 
ruKiT K, lOOGeV. This translates into that the scalar contribution to the fermionic self-energy is not sensitive to the 
infrared behavior and does not require hard-thermal loop resummation. 

By casting our studies in terms of the real time, in-medium effective Dirac equation, we were able to study the real 
time evolution of wave packets of collective excitations and extract unambiguously the mean free paths or attenuation 
lengths for the different collective excitations. This approach transcends any approximation that relics on a diffusive 
description of transport. 



Our main results are summarized by the expressions: (3.36) for the scalar contribution to the damping rate and 
(4.9) for the mean free paths for the soft collective excitations with k ^ g^T and away from the plasmino minimum 
k K, 0.4 Meff- The main observation is that the difference of mean free paths for the fermion and plasm ino branches 
for k ^ g^T and away from the minimum of the plasmino branch is approximately given by Eq. (4.10), i.e., mainly 



determined by the different scalar contributions to the damping rate of the collective excitations. The uncertainties 
near the region of the plasmino minimum are related to the infrared sensitivity of the fermionic damping rates and 
require a careful analysis of the infrared region of the HTL resummed self-energies, such a study is currently under 
way. 

There is potentially an important exception to the validity and robustness of our results in the case of extensions of 
the Standard Model that include scalars that are too light. Light scalars are common in supersymmetric extensions 
and extensions with more than one doublet. They are welcome because they tend to increase the strength of the first 
order phase transition [^. If light scalars are present, their contribution to the fermionic self-energy will have to be 
understood in detail and assess in each particular model whether an HTL resummation of the scalar propagator as 
well as that of the fermion in the loop is necessary. Clearly such cases must be analyzed in detail for the particular 
models. 

A very important limitation of our results is the knowledge of the damping rate of the plasmino near the region 
where its group velocity vanishes. As recognized by the authors p^- [2^j3^ of the original work, in this region there 
are further contributions to the plasmino damping rate that must be understood in detail. 

We think that the novel results found in this article for the transport of chirality by soft collective excitations reveal 
new features which are distinct from those previously obtained for hard excitations and could have an important impact 
in mechanisms of local baryogenesis. Whereas the main focus of this article has been a quantitative understanding of 
chirality transport by soft collective excitations and a detailed understanding of the gauge and scalar contributions 
to the transport coefficients, obviously the next step is to provide a quantitative assessment of the impact of these 
new results on the baryon asymmetry within particular extensions of the Standard Model. This study is currently 
underway. 
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FIG. 1. Dispersion relations: the upper correponds to the fermion branch and the lower to the plasmino branch. The dashed 
hne is the light cone. 



0.08 



0.06 



3~ 



u 0.04 



0.02 




FIG. 2. rj_/Me// vs k/M^ff for g — 0.3. The solid line corresponds to the fermion branch and the dashed line to the 
plasmino branch. 
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FIG. 3. 32TTr'^{k)/{y^Meff) vs k/M^ff for m — T. The soli d line is for the fermion branch and the dashed line for the 



plasmino branch. The scalar contribution rj;(fc) is given by Eq. (3.36) in the text 
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FIG. 4. [A+(fe) - \-{k)\Me.ff vs k/M^jf for g = 0.3, y = lO""* and m = T. The peak near k « 0.4 Me// is a consequence 
of the vanishing of the plasmino group velocity, the region near this peak is not trustworthy as the expression for the gauge 
contribution to the dam ping rate of the plasmino is not valid in the region where its group velocity vanishes as explained in 
the text below Eq. ([j.26|) . 
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